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We investigate adiabatic charge pumping in disordered system in one dimension with open and 
closed boundary conditions. In contrast to the Thouless charge pumping, the system has no gap 
even though all the states are localized, i.e., strong localization. Charge pumping can be achieved 
by making a loop adiabatically in the two-dimensional parameter space of the Hamiltonian. It is 
because there are many function- like fluxes distributed over the parameter space with random 
strength, in sharp contrast to the single ^-function in the pure case. This provides a new and more 
efficient way of charge pumping and polarization. 



INTRODUCTION 

The quantum dynamics of the charge in insulators is 
a rich and non-trivial issue even in the d.c. limit. Of 
course, the d.c. conductivity vanishes, but this does not 
mean that charge motion is frozen in the insulators. One 
example is the ferroelectricity and electric polarization 
in insulators. The classical definition of the electric po- 
larization P = J drrp{r) ( p(r): charge density) fails 
for the extended Bloch wave functions since the f is un- 
bounded. This difficulty is avoided by considering in- 
stead the new definition for the difference of the polar- 
ization 0,0, 0,0,11, 



AP 



T 1 dP 
dr—, 
dr 



(i) 



where the change of the polarization between the initial 
and final states are given by the integral of the polariza- 
tion current during the adiabatic change of the param- 
eters Q = (Qi, Q2, • • • , Qn) such as the atomic displace- 
ments. One can usually choose the initial state with the 
inversion symmetry without electric polarization, and (pQ) 
uniquely determines the polarization of the final state 
of our interest. Here by using j£ = ^j-^S the /> 
component of AP is expressed as 



Jc dQ 



(2) 



with the path C specified by Q = Q(r). What is found 
by Rest a [3] and King- Smith and Vanderbuilt 0] is that 

dP 

—j^r can be represented by the Berry phase, which fits to 
the first-principle band calculation. Then the question 
arises, "Is there any path (C) dependence of the polar- 
ization ?". It is clear that there is no parametrization 



dependence since we consider the adiabatic change, but 
the different path C\ and C2 in the Q-space might lead to 
the different values of AP. This is related to the single- 
valueness and Chern number of the Bloch wavefunction. 
Onoda et at. addressed this problem as follows by analyz- 
ing the one-dimensional two-band models characterized 
by the three dimensional Q 0,0- There appears a singu- 
lar line, i.e., string, in the Q-space corresponding to the 
trajectory of the monopole (band crossing point) as the 
momentum k moves, which acts as the "current circuit" 
to produce the "magnetic field" via the Biot-Savart 
law. Away from the string, the system is always gapped, 
and insulating. When the adiabatic change of the param- 
eter Q makes a loop enclosing the string, the charge is 
pumped during this process, which is quantized to be an 
integer multiple of e since the strength of the "current" 
is quantized. This is a realization of the quantum charge 
pumping first proposed by Thouless [8[. This analogy 
to the magnetostatics says that the polarization is path 
independent as long as the loop C = C\ + (— C2) does 
not enclose the string. This is usually the case because 
the change of the parameter Q is rather small and we 
need the huge variation of Q to enclose the string, i.e., 
gapless states. In other words, the ferroelectricity can be 
regarded as "a fraction of the quantum charge pumping" . 

This charge pumping can be regarded as the rigid shift 
of the wavefunction due to the change in the external 
parameters such as the atomic positions. It is natural 
when the wavefunction consists of the Bloch states ex- 
tended over the whole sample, and then the quantum 
interference pattern is modified by the external parame- 
ters. One can estimate roughly how much the charge is 
pumped as below. Let the dimension of the parameters 
be the energy. Then the distance between the physically 
realized set of parameters and that of gapless states, i.e., 
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string, is the energy gap Eq- When A be the change of 
the the parameters in unit of energy, the angle subtended 
by this segment in parameter space is roughly estimated 
as 2ttA/Eg- Since the 2ir winding corresponds to unit 
charge e shift by one lattice constant a, i.e., Po = ea, the 
polarization is roughly given by P = eaA/Ec [6[. There- 
fore one can enhance the dielectric response by reducing 
Eg, or enlarging A. One possible method to reduce Eq 
is by introducing disorder, by which even the gapless 
insulator can be realized. However, the electron wave- 
functions are no longer the extended Bloch states in the 
presence of the disorder, and one needs to worry about 
the localization, i.e, Anderson localization. When all the 
states are strongly localized, one can not transmit the 
phase information though the sample, and hence can not 
expect the charge pumping either. On the other hand, 
as can be seen from the above explanation, the charge 
pumping is closely related to the topological nature of 
the wavefunctions in the parameter space, which is ro- 
bust against the disorder to some degree. For example, in 
the two dimensional electron systems under strong mag- 
netic field, there occur discrete extended states •protected 
by the topology. Namely, the Chern number is carried by 
the extended states only, which is not destroyed by the 
weak disorder. Therefore similar problem arises even for 
the one-dimensional systems which we study below. Niu 
and Thouless [9] was the first to study the stability of the 
charge pumping against the weak disorder by a topolog- 
ical argument. Even though all the states are localized, 
the charge pumping was shown to be unchanged as long 
as the gap remains finite. However, questions still remain 
such as what is the physical mechanism of the charge 
pumping through the localized states and what happens 
when the gap collapses. 

In this paper, we study the effect of the disorder on 
the charge pumping and dielectric response in a one- 
dimensional model for insulators. Combining the numeri- 
cal simulations and analytic considerations, we reveal the 
physical picture of the charge pumping by the localized 
states both in the case of open boundary condition and 
periodic/twisted boundary condition. The former one is 
more relevant to the experimental situation such as the 
FeRAM (ferroelectric random access memory) where the 
leads are attached to the thin film of insulators, while 
the latter is more appropriate to see the role of topol- 
ogy. We have published a Letter focusing on the open 
boundary condition [10[. This paper provides the full 
details of the formulation, calculations, and additional 
results on the open system, as well as new results on the 
periodic/twisted boundary condition. 

The plan of this paper follows. In section II, a model 
for disordered insulator showing the charge pumping is 
introduced. Its analysis with the open boundary condi- 
tion is given in section III, including the detailed descrip- 
tion of the resonant tunnelling mechanism. In section IV 
given the analysis of the model in the three dimensional 



parameter space including the angle a for the twisted 
boundary condition, and the role of magnetic monopoles 
in this space. Section V is devoted to the conclusions. 

A MODEL FOR DISORDERED INSULATOR 

The minimal model for ferroelectrics is given by the 
following ionic dimer model; 

1 L 

#pure = - - ( tn - n - ~ (-YQl) ( C Ul°i + h ' C ') 

2=1 

L 

+ J2(-YQi4a- (3) 

2=1 

Here, ci and c\ are the annihilation and creation opera- 
tors of the electron at the site i = 1, • • • , L with L being 
the number of sites in the system. For open system at- 
tached to leads at the both ends, c^+i and c^ +1 repre- 
sents those operators in one of the leads, while for closed 
system, they are understood as c\ and c\. t n . n . is the 
transfer integral. Qi and Q2 represent the alternations 
of the local ionic level and the bond dimerization, re- 
spectively. The spin degree of freedom is omitted for 
simplicity. Then, of our interest is the half-filling case 
relevant to the ferroelectrics. Although this model might 
look special, it represents the two essential feature of the 
ferroelectrics, i.e., (i) the two species of the ions charac- 
terized by the level alternation Qi and (ii) relative shift of 
the atomic positions described by the dimerization Q2. It 
corresponds to the ferroelectricity in e.g. BaTi03, where 
Ti and O are dimerized to produce the polarization . 
It can be also applied to the quasi-one dimensional fer- 
roelectric materials such as organic charge transfer com- 
pounds TTF-CA [12] and (TMTTF) 2 X [13[. 

The Hamiltonian H pure under the periodic boundary 
condition yields two bands [6] 

e±(k) = ±^tl n .cos 2 k + Ql + Q 2 2 sm 2 k. (4) 

Experimentally, the parameters Q = (Qi,Q2) can be 
controlled by applying the electric field E along the polar- 
ization direction and the pressure p as follows. Electrons 
at high and low deinsity sites shift in relatively opposite 
directions, as shown in open arrows of Fig. [1] (b), chang- 
ing the dimerization Q2 by SQ2 oc eEQ\. Simultane- 
ously, within each dimer, a level difference Q\ changes by 
SQi oc —eEQ2 as illustrated by black arrows in Fig.[T](b). 
Therefore, the electric field E mainly controls the angle 

= arctan(Q 2 /Qi). (5) 

Applying the pressure, one can increase the hybridiza- 
tion, and then reduces the ratio Q/t n .n. with 

Q = ^Qi + Qi (6) 
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To discuss the quantum relaxor behavior, we introduce 
the on-site random potential V{ to the Hamiltonian given 
by ©; 



H = H, 



pure 



E ViC i Ci 



(7) 



The type of the random distribution takes either a uni- 
form distribution or an alloy model, which are shown 
in Fig. [21 In the following, we study effects of on-site 
disorder on dielectric properties in both open and closed 
systems, particularly focusing on the topological aspects. 
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FIG. 1: (Color online) Ionic dimer system sandwitched by the 
leads. 



(a) 



Probability 



(b) 



5v 



Probability 



■v 



) 

V v : 



8v 



-V V V; 



FIG. 2: (a) Uniform distribution of the on-site random poten- 
tial, (b) Alloy model for the on-site random potential, which 
mimics effects of substitution in alloys. 



OPEN BOUNDARY CONDITION 

In this section, dielectric property of disordered insula- 
tors connected to the leads is studied theoretically. This 
corresponds to the realistic situation in FeRAM devices, 
where the polarization is measured by the integrated cur- 
rent flowing in the leads. Especially in nano-scale sam- 
ples of the FeRAM, the quantum nature of the polariza- 
tion is expected to play essential roles via novel quan- 
tum interference as discussed in the introduction 0, @] . 
Namely, a theory for the nano-scale capacitance is called 
for where the quantum coherence is attained all through 
the sample. ^From the viewpoint of the applications, it 
is required to achieve (i) a magnitude of the polarization 
larger than 10 /iC/cm 2 , (ii) a smaller leak current than 
0.1-1 /iA/cm 2 , and (hi) a dielectric constant larger than 
300 [141 ]. By introducing disorder, one can also suppress 
the dissipation by leak current due to the strong local- 
ization effect [3, [l5[ in addition to the reduction of the 
gap and enhanced dielectric response. In fact, a similar 
disorder-driven enhanced resp onse has been found in re- 
laxor ferroelectrics [lg, LL7|, LL8| and pinned charge density 
waves [HI, although its mechanism is due to a classical 
mesoscopic cluster formation analogous to spin glass sys- 
tems. 

In the disordered case with all the states being local- 
ized, the charge transfer occurs through the resonant tun- 
nelling. The quantum nature is topologically protected. 
Crossover from the quantized charge pumping in bulk in- 
sulator [8[ to that of quantum dot [20] is clarified in the 
control parameter plane as a function of sample size and 
disorder strength 2l|, [H, Q . Then, we theoretically pro- 
pose the enhancement of the dielectric response of insula- 
tors by disorder in nanoscopic/mesoscopic multichannel 
systems, with the required time for the dissipationless 
adiabatic charge displacement being estimated. Appli- 
cation of this phenomenon to memory devices, like fer- 
roelectric random access memory (FeRAM) [14[, is also 
discussed. 



Single-channel problem for uniformly distributed 
random potential 

We consider an insulating electronic system sand- 
wiched by leads (electrodes) as shown in Fig. 2(a). For 
simplicity, we take a one-dimensional (single-channel) 
model, but the extension to higher dimensional (multi- 
channel) cases is straightforward as discussed later. 

We take the total Hamiltonian 



Htot — H + i^lead 
with H and #i ea d being given by Q and 

\i=0 i=Lj 



(8) 



ffiead = E + E + k - C ) ' ( 9 ) 
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respectively. This model is schematically shown in 
Fig. [I] (b). The Green's function Gi^(e) for the above 
model is readily obtained from the recursion formula in 
the form of the continuued fraction [2J]. Here, we can 
concentrate on the case where the chemical potential is 
located at the zero energy. 

We adopt the Landauer-Biittiker formalism 



24, m 



where the sample is regarded as a scatterer characterized 
by the scattering S'-matrix 



S 



r t 
t* r' 



(10) 



Here, the reflection coefficients r and r', and the trans- 
mission coefficient t (see Fig. [1]) can be calculated from 
the Green's functions [24]. The transmittance through 
the sample and the reflectance at the both ends are ex- 
pressed as 



T 
R 



>|2 



1-T, 



(11) 
(12) 



respectively. Then, we employ the Brouwer's formula [26L 
I27I ]. which has been successfully applied to the charge 
pumping in metallic quantum dot systems [20]: the 
charge Aq pumped from the left to the right during an 
adiabatic change of parameters Q along a path C is given 
by M 



Aq = ej^-lm (r* V Q r + tV Q t*) . 



(13) 



Therefore, even in the perfectly reflecting case t = 0, \r\ = 
1, charge can be pumped by controlling the phase cp of 
r = 0. 

Let us start with the pure case (vi = 0). Without the 
random potential v^, the bulk system has an energy gap 



E G0 = 2jQj + Mm{tl^Q 2 2 }. 



(14) 



This gap closes at Q = , where the metallic conduction 
occurs. This yields a vortex center of the reflective co- 
efficients r and r' . By means of the continuued- fraction 
expansion of the Green's function, an analytic form of r 
in the thermodynamic limit L — » 00 is obtained as 



Q\ + tn.n.Q2 + Qy^.n. + gj + 2zt n . n .Qi 
Qi + *n.n.Q 2 + QVtln.+Ql " 2z* n .n.Ql 



where Q = vQi+Ol- The panels (al) and (bl) of 
Fig. [3] show the reflectivity R = \r\ 2 for L = 10001 and 
the phase cp of r for L — > 00 without randomness. Here 
the phase winds by 2tt around Q = (0, 0), forming a "vor- 
tex" at which r = and \t\ = 1. In the thermodynamic 
limit L — » 00, R = 1 and T = hold except at the vortex 
center. For finite-size systems, the transmittance T be- 
haves as e~ L ^° with £ = tn.n./^GO- Therefore the size 



of the region in the Q space with a large T is of the order 
of t n . n jL for a large L. Namely the vortex corresponds to 
the gapless case, where the extended state at the Fermi 
energy carries the charge and causes the perfect trans- 
mittance. Remarkably, this perfect-transmission point 
Q = governs topological properties of the system in 
the whole Q plane in a non-local way. When we adiabat- 
ically change the parameter Q along a cycle around Q = 
at which the transmittance T vanishes, the charge e is 
pumped and the polarization changes by ±2ea (a: lat- 
tice constant) according to the Brouwer's formula (p~3|) . 
The vortex is almost isotropic at least in the vicinity of 
its core. Then, a pumped charge due to a small change 
A of Q is expressed as q ~ (</>/27r)e ~ (|A|/i^)e, where 
4> represents a change in the polar angle of the vector Q. 
Therefore one can enhance the pumped/displaced charge 
by enlarging the angle <\> subtended by A around the gap- 
less point. This quantized charge pumping [8] through 
the adiabatic cyclic change of Q is consistent with the re- 
sults previously found in the periodic system [6] by using 
the Berry-curvature formulation of the electric polariza- 
tion m. 

Now we turn on the disorder in order to reduce the 
gap Eq- In the Q plane, the region of \Q\ < v be- 
comes gapless, while the Anderson localization seriously 
affects the transport properties [15|, S] . In one dimen- 
sion, the effect is pronounced and all the states are local- 

U El, 



ized 



On the other hand, the total vorticity 
N v is an integer topological number, and is robust under 
a continuous change of parameters including the disor- 
der strength. Therefore, even with disorder, N v remains 
unity. Namely, there should exist at least one vortex 
(r = 0) in the Q plane; thus perfect transmittance \t\ = 1 
occurs even though all the states are strongly localized. 
This is a remarkable consequence from the quantum and 
topological property of the system, and is in sharp con- 
trast to usual classical tunnelling through disordered in- 
sulators. 

In Figs. El (a2) and (b2), we show our numerical results 
of R and cp for L = 101 with a uniform random distribu- 
tion Vi G [— v,v] for the disorder strength v = 0.25t n .n.- 
The vortex center of the perfect transmittance T = 1 and 
R = shifts in the Q space. Besides, an anisotropy devel- 
ops in the shape of the region of relatively high transmit- 
tance T. We further calculate R and cp for larger systems, 
which are shown in (a3) and (b3) for L = 201 and in (a4) 
and (b4) for L = 401. There, it is evident that the vor- 
tex core, which is almost isotropic without the disorder, 
rapidly evolves into the highly anisotropic one. This is 
associated with an increase of the ratio L/^by increasing 
the system size. Figure |4] shows the local density of states 
Di of site i{= 1, • • • , L), which has been calculated when 
Q is located at the vortex center Q c corresponding to 
each case of L = 101, 201 and 401 for v/t = 0.25. When 
L is equal to 101 or smaller, the state at this energy is 
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FIG. 3: (Color) Reflectance R = |r| 2 , and the phase ip = 
arg r. (al), (bl): clean bulk system. (a2), (b2): disordered 
system of finite size L — 401 with the random potential of the 
strength v/t n .n = 0.5. The color code given in (al)/(bl) also 
applies to (a2,3)/(b2,3). In the white region in (al,2,3), R is 
the unity within the accuracy of 10 -4 . 



extended over the sample. However, with increasing L, 
the state is almost localized in the middle of the sample. 
The spatial extension of the wave function, i.e., the lo- 
calization length £, can be explicitly evaluated from the 
second moment of Di as the inverse participation ratio; 



6pr = (E^) 2 /E^ 2 - 



(16) 



We obtain £ IPR = 51.4, 80.8, 76.4 and 84.7 for L = 101, 
201, 401 and 501, respectively around the vortex center, 
indicating that £ipr almost saturates about 80 sites for 
v/t 0.25. 

In this anisotropic "wing", ip changes rapidly. The 
width Qw of the anisotropic wing decays exponentially 
as exp(— L/£) with increasing L, as shown in Fig. [5j 

In the rest of this subsection, we comment on the re- 
lation of the charge displacement discussed above to the 
charge pumping through a metallic dot [20[ weakly con- 
nected with two leads. In the metallic dot systems, a 
charge is pumped from one lead to the other by chang- 
ing the heights of potential barriers between the dot and 



(a) 



(b) 




(c) 



FIG. 4: Local probability Di of the state at the vortex for 
(a) L = 101, (b) L = 201, and (c) L = 401 in the presence 
of uniformly distributed random potential v% £ [—v, v] with 

v — 0.25£ n .n. . 



the neighboring leads [26[. Since the system is metal- 
lic, the pumped charge is not quantized in this case. 
Here, the main source of the pumping is a finite trans- 
mitt ance but not the phase of reflection coefficient. With 
increasing the system size, it crosses over to the opposite 
regime with vanishing transmittance studied in this pa- 
per. In the present theory, the charge pumping comes 
from the vortex of the reflection coefficient. Therefore, 
it is topologically protected against disorder. This topo- 
logical constraint guar antees the applicability of the res- 
onance tunnelling [33], only near the vortex. 
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100 200 300 400 

L 

FIG. 5: The minimum width SQ for R determined at R = 0.3 
as a function of L. 



Resonance tunnelling 

In order to understand transport properties of such 
one-dimensional disordered systems, it is sometimes use- 
ful to consider an effective model with the potential hav- 
ing high double peaks. In this model, well-defined local- 
ized eigenstates exist between the two potential peaks, 
and any transport between two ends of the system oc- 
curs via the localized states through tunnelling. While 
such tunnelling has exponentially small probability, it oc- 
curs when the Fermi energy of the leads is equal to one 
of the eigenenergies of the localized states, namely when 
the resonance takes place. This picture of transport in 
disordered systems is called "resonance tunnelling" [331 ]. 

Here we describe the resonance-tunnelling theory, fol- 
lowing Ref. |23, and will see whether it fits with the 
present model with disorder. We define an effective 
model of "resonance tunnelling" as described by the 
Schrodinger equation ^"{x) + (k 2 — V(x))^(x) = where 
h 2 k 2 /2m is the particle energy, and the effective poten- 
tial V(x) has two peaks as shown in Fig. [6l It is not 
trivial whether the numerical results of our model ([3]) fits 
well with this resonance-tunnelling picture. The present 
model has many valley and peaks, due to on-site random- 
ness, and is not similar with Fig. [6l Nevertheless, when 
the localization length £ is much shorter than the system 
size, it is indeed the case as we see by fitting our numer- 
ical results well by the picture of resonance tunnelling. 
The basic reason for the applicability of resonance tun- 
nelling is that the topology guarantees an existence of a 
perfect-transmittance point, i.e., the resonant tunnelling 
in the two-dimensional parameter space (Qi, Q 2 ), in spite 
of the complexity of the potential shape in the model. 

First let us consider an open system attached to two 
ideal leads which are connected to reservoirs. The Fermi 
energy E F is tuned at one's disposal, and we assume that 



X1 X'1 X"1 X2 X'2 X"2 X3 

FIG. 6: Potential profile for the effective model of resonance 
tunnelling. Heights of the two potential peaks are character- 
ized by Si and $2. 

the potential peaks are much higher than Ep. An effec- 
tive wave number k(x) = \Jk 2 — V(x) becomes imagi- 
nary inside the potential peaks, and the particle tunnels 
through the peaks. One of the eigenstates at the left lead 
(x < x\) is given by ^\{x) ~ exp(z/ci(x — #1)), where 
k\ = k(xi). The time-reversal symmetry of the system 
requires that the complex conjugate S&i*(x) is also an 
eigenstate. At the potential valley around x = X2 we can 
also consider eigenstates with ^2^) ^ exp(z/c2(x — #2)) 
and S&2*(%)i where ki = k(x2). Thus general eigenstates 
around x = x\ and around x = X2 are given as ty(x) = 
Ai^iO) + and V(x) = A 2 ^ 2 (x) + £ 2 #2*(a), 

respectively. Next we introduce the transfer matrix ©1 
from X2 to x\ 




where Bi is a 2 x 2 matrix. Time-reversal symmetry 
requires the transfer matrix to be of the form 

a - (h\ V cosh(Si)e tol sinh(Si)e^ \ 

UJ V sinh (Si)e" i/31 cosh(Si)e- iai )• [ ^ } 

In a semiclassical theory we have 

ai = — / k(x)dx — / k(x)dx, (19) 

(3 1 = / k(x)dx + / k(x)dx. (20) 

Si characterizes the height of the potential, as the trans- 
mission coefficient t\ is given by t\ = e m sechS'i. Because 
we assume that the peak is high, which means Si ^> 1. 
The transmission probability |ti| 2 is then approximated 
as |ti| 2 - 4e~ 2Sl < 1. 

Similarly, we introduce the transfer matrix for the sec- 
ond peak as 
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and we have 



e 2 = ^ 



cosh(S 2 )e ia2 sinh(S 2 )e i/32 
sinh(S , 2 )e-^ 2 cosh^e"*" 2 



• (22) 



The phases a 2 and /3 2 can be written similarly to a± 
and ySi, where S2 ^> 1. We assume that Si, ol{, and 
(z = 1, 2) are smooth functions of the parameters Qi, Q2 
and Ep. A transfer matrix for the entire system is 
given by 



in = e 



11 u 21 
21 C/n 

i(«i+« 2 ) ( cosn Si C oshS 2 + sinh Si sinhS 2 e 



(23) 



i2i = e 



(sinh Si cosh S2 + cosh Si sinh S26 



(24) 



(25) 



r x 2 

J x'J 



uj = pi - f3 2 - oil - a 2 = 2 / k(x)dx. (26) 



/From unitarity, it follows that |0n| 2 — \0 2 i\ 2 = 1. We 
assume that the Fermi energies of the two leads are iden- 
tical, which means ks = fci- For a plane wave incident 
from the left lead, let r and t denote the reflection and 
transmission coefficients, respectively. Similarly, for a 
plane wave from the right lead, we define r' and t' as 
well. We then obtain 



^21 



/0 1U t = t' = I/flu, r = -02*1/011. (27) 



It implies the unitarity |r| 2 = |r / | 2 = l — \t\ 2 . It also 
satisfies t = t' and rt* + r'*t = 0, as is required from 
time-reversal symmetry. 

We now apply this framework to fit our numerical re- 
sults. Because r is written as 



r = e 



•A tanh Si + tanh S 2 e l{ 



(28) 



1 + tanh Si tanh5 2 e iw ' 
a condition for a total transmission, r = 0, is given by 

Si=S 2 , (29) 

f X 2 

u = 2 k(x)dx = (2n + 1)tt, (30) 



where n is an integer. The latter condition is equiva- 
lent to the Bohr quantization condition, that the state 
localized around x 2 be an eigenstate with its eigenenergy 
equal to Ep. In other words the total transmission oc- 
curs at resonance. For fixed Ep, the two conditions (|29|) 
and ([3Q|) define isolated points in the Qi-Q 2 plane. Let 
P denote one of such points of total transmission: r = 0. 
One can easily see that around the point P, the phase of 
r and that of r' wind by ±2tt, as is schematically shown 
in Fig. EKal) (a2). As is seen from the Brouwer formula 
(fT3]h the phase windings of r and r' correspond to (2tt 



times) the amount of charge pumped into the system 
through the left and the right ends, respectively. Thus, 
by going around the point P, one unit charge is pumped 
from the left lead to the right. For illustration, let us 
consider a clockwise cycle around the point P in Fig. [71 
(al)(a2). This pumping from left to right is analogous 
to the "bicycle pump" [28]. The two potential peaks 
correspond to two gates to control the pumping. If the 
system crosses the line uj = (2n + 1)tt at the Si < S 2 side, 
r undergoes 2tt phase changes. It means that tunneling 
occurs through the left potential peak due to resonance, 
corresponding to the opening of a "left gate" and a unit 
charge flows in. After that the left gate closes, and the 
right gate opens in turn, when the system crosses the line 
uj = (2n + 1)tt at the Si < S 2 side, and r' undergoes — 2tt 
phase change. The charge is pumped to the right lead 
after one cycle. Thus the overall motion of the charge is 
as shown in Fig. [71 (a3). For the pumped charge to be 
quantized, the cyclic process should be sufficiently slow 
to be regarded as "adiabatic". Otherwise the particle 
cannot tunnel through the potential barriers. Below let 
us make this physical picture more explicit in relation to 
our numerical results. 

In our numerical results in Fig. [3jb2)-(b4), the phase- 
winding point P corresponds to the perfect transmission. 
When we go from Fig.[3](b2) to (b4), L/£ increases. In the 
resonance- tunnelling picture, it means that Si increases 
and that the "resonance- tunnelling wing" narrows, as is 
seen by comparing Fig. 0[al) and (bl). It is exactly 
observed in our numerical results. 

For large L/£, the transmission probability \t\ 2 is typ- 
ically ~ e -2(Si+s 2 ) ^ e -2£/£. Because we have Si ~ S 2 
near the resonance, we estimate Si ~ S 2 ~ (L/2£). In 
that case, the change of the phases of r and r' occurs 
abruptly at around \uj — (2n + 1)tt | ~ e~ L ^. Along the 
line uj = (2n + 1)tt, |r| is given by \r\ = tanh (Si — S2), 
which becomes appreciable only when |Si — S2I °° 1- This 
is found only in a small region A/t n . n . ~ £>/L. When 
L/£ ^> 1, because the tunnelling rate Qw is exponen- 
tially small (Qw ~ t n . n .e~ L ^), charge pumping requires 
an exponentially long time r ~ (^/^ n .n.)e L// ^. This is 
a time scale which gives a criterion for adiabatic charge 
pumping in this system. If one changes the parameters 
faster than this time scale, the pumped charge is reduced 
with an exponential factor e - ^/^ as a nonadiabatic cor- 
rection, where a; is a frequency of the change the external 
parameters [34j . 

Generally, there occurs no perfect transmittance point 
by tuning only two parameters since the effective two 
barrier model for the resonant tunnelling applies only to 
a limited region of the random system and not through 
a whole sample. In the present case it is guaranteed by 
the topological constraint that the phase of r (r') should 
wind by (— )2tt for a large cycle far from Q = 0, well 
within the gapped region [9[. When \Q\ is larger than the 
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(a1 ) phase of r (small 81,82) (a2) phase of r ' (small &,&) 




(a3) charge motion 

[A-*- B CJ [C-^D-^A] 




(b1) phase of r ' (large Si,S2) (b2) phase of r 1 (large Si,S2) 




FIG. 7: Schematic contour mapping of the phases of the re- 
flection coefficients r and r' in the Q1-Q2 plane. (al)(a2) 
correspond to small Si, &, and (bl)(b2) to large Si, S2. The 
phases of r and r' wind by 2tt during the cycle around the 
"vortex" P. When we go along the arrows indicated in the 
figures (al)(a2), the phase of r (r) increase by 2ty (—2tt). (a3) 
illustrates the motion of the charge in the cycle. Note that 
the 2tt phase change of r (r ; ) is associated with pumping one 
unit charge into the system through the left (right) end of the 
system. 



Multi-channel problem in alloy model for on-site 
random potential 

In the rest of this section, we consider a stronger ran- 
domness by employing the alloy model (Fig. [2] (b)) with 
Vi = Si(v + Svi) with Si = ±1 being a random sign and 
Svi/v G [—0.025,0.025] a uniform random distribution, 
instead of the uniformly distributed random potential 
(Fig. El (a)). 

Let us start with the single-channel case. Then, for 
a fixed disorder strength v, the shape of the vortex core 
of the reflection coefficient rapidly evolves from isotropic 
to anisotropic with increasing the system size L, as in 
the case of the uniformly distributed random potential 
discussed in Sec . This tendency is generic and also real- 
ized for a fixed system size L with increasing the disorder 
strength v. In Figs. [H (a) and (b), we show the reflectance 
R and the phase tp in the Q space for v/t n , n , = 1.0 and 
a small sample size L = 25. Here, the region in the Q 
plane where the gap collapses and vortices distribute ex- 
pands because of the stronger disorder and thus the much 
shorter localization length £ ~ 3 or 4. The transmittance 
T is typically obtained as ~ 10 -9 , which is practically 
negligible except at the vortex core. 



(a) R (b) <p/2n 




QiAn.n. QlAn.n. 



energy scale of the disorder, a finite gap (~ |Q|) opens. 
The wing will then become as wide as Qw ~ \Q\- Cor- 
respondingly, the typical time scale r ~ Tt/Qw ~ h/\Q\ 
becomes smaller, and the adiabaticity condition is easily 
satisfied. In such case, however, the dielectric response is 
not enhanced. The localization length becomes as long 
as the system size, and the pumping is accomplished 
through extended states, not by resonance tunnelling. 
Remarkably, this charge pumping in the gapped region 
is governed by the vortex which is located deep in the 
disordered regime (£ <C L). In other words, the charge 
pumping in the gapped regime (£ ~ L) is smoothly con- 
nected to that via resonance tunnelling in the disordered 
regime < L). We note that S\ and £2 are regarded as 
effective parameters, although the real potential is much 
more complicated than the two-barrier structure. 



FIG. 8: (a) Reflectance and (b) the phase ip of r in disordered 
system with L = 25 and v/t n . n . = 1.0 for the alloy model with 
the on-site random potential given by Fig. [2] 

Now we consider a system that consists of many chan- 
nels, each of which is described by the present disordered 
alloy model but with a different profile of random poten- 
tials. Such configuration can be realized in thin films 
of ferroelectrics. Then, we can design the pattern of 
the phase cp in the Q plane by tuning the disorder. For 
the alloy model with the random on-site potentials, vor- 
tices are mostly located around Q\ = ±v/2 with the 
"wing" almost along the Q2 axis. Therefore, we can en- 
hance the dielectric response, if we can tune the disor- 
der strength and choose a sample where the Q point of 
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the system is located inside the "wing". In particular, 
when the transmittance T is negligibly small, the en- 
hancement factor of the dielectric response is given by 
(dip/ 89) disorder/ (dip/ dO) pure from ([T3]h since the electric 
field is proportional to — arctan((5i/Q2)- We calcu- 
late this enhancement factor for this multichannel sys- 
tem. Increasing the number N of channels. Figures [9^a) 
and (b) represent the Q dependence of the enhancement 
factor in the case of v/t n , n , = 1.6 and 2.4, respectively, 
for L = 25 with N = 102. The main structure in this 
map is almost saturated up to N = 10 2 . These results 
reveal that around Qi ~ ±^/2, the dielectric response 
is significantly enhanced by a factor 30 ~ 40 compared 
with the pure case. Even for the thin film with a square 
shape of a linear dimension larger than 50A, which cor- 
responds to N = 10 2 , the disorder-induced enhancement 
of the charge transfer rate should be robust. Then, the 
applied electric field necessary for switching the polar- 
ization is reduced by this enhancement factor. If we re- 
quire the response time r ~ e L ^/t n . n . of the order of 
10 -9 s, we obtain e L ^ < 10 6 with the assumption of 
t n .n. ~ 10 15 s _1 . This also assures a negligibly small 
transmittance T ~ e~ L ^ ~ 10 -6 , and thus the small 
leak current and low dissipation. 




Ql An.n. Qi An. n. 



PERIODIC/TWISTED BOUNDARY CONDITION 

The adiabatic charge pumping in the presence of the 
substrate disorder was considered by Niu et al.[9}. They 
showed that the adiabatic charge transport is still quan- 
tized as long as the excitation gap between the high- 
est occupied state (HOMO) and lowest unoccupied state 
(LUMO) does not vanish due to the substrate disorder 
and the many-body interaction in the thermodynamic 
limit. In the case of open boundary condition, it is easy 
to imagine the meaning of the charge pumping, namely 
the charge transport from one end to the other. In the 
case of periodic boundary condition, the charge pumping 
after one-cycle means that the electronic wavefunction 
shifts by the minimum number of lattice sites so that the 
final wavefunction is the same as the initial one. A car- 
toon picture of this process is given by Fig. [10l We note 
that in the absence of the disorder, the adiabatic charge 
transport has been related to the field-theoretical model 
of the one-dimensional chiral anomaly [37j. 




FIG. 10: A cartoon picture of the charge transport on the 
one- dimensional ring. 



Formalism 



FIG. 9: Relative dielectric response in the presence of the 
disorder compared with the pure case. The disorder strength 
is v/t n . n . = 0.8 for (a) and 1.2 for (b). Averages are taken 
over 10 2 random disorder configurations. Inside white bands, 
there occurs a gradual sign change in the dielectric response. 



Possible experimental realization of this quantum- 
mechanical disorder-induced enhancement of the dielec- 
tric response, namely the quantum relaxor, has also been 
proposed [10] for thin films of solid solution systems like 
Pb(Feo.5Nbo.5)03 and Pb (Sep. 5 Nbn 5 ) O3 prepared with 
an adequate slow-anneal process [35|, [36[. 



We consider a one-dimensional charge-density-wave 
system with disorder, as described in (3) and (7). In 
contrast to the open system with two leads as considered 
in Section III, we deal with a closed system in a ring 
geometry. Instead of a periodic boundary condition as 
is usually employed for the ring, we introduce a twisted 
boundary condition; we impose that the phase gained 
by running around the ring once is e ia . This twisted 
boundary condition is realized by applying a magnetic 
flux through the one-dimensional ring. 

The reason for introducing the twisted boundary con- 
dition is to reveal the topological nature of the charge 
transport in the system. In the pure case, the lattice 
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momentum k is a good quantum number, and the topo- 
logical nature is manifest in the three-dimensional space 
as monopoles and antimonopoles. However, 
k is no longer a good quantum number in the disordered 
case, and instead the flux a plays a similar role to k. 
Therefore, instead of (Qi,Q2,&), we will be interested 
in the adiabatic process such that the system is changed 
due to the slow variations of the parameters (Qi, Q2, a), 
which spans a 3-dimensional parameter space with the 
property a G [0,27r]. The Hamiltonian with the flux a 
can be rewritten as 



H(a) 



^n.n. 



N 



^(e %a/N c]c J+1 ^h.c.) 

N 

Q\ E(- 1 )' c J c j 

3=1 

Q N 



^1 Y^{-iy{e^ N c)c 3+1 + h.c.) + H31) 



is the uniformly 



where c/v+i = c\ and V 
distributed disorder potential. 

In the tight binding model, the polarization operator 
P has the following form 



(32) 



where Rj is the position at site j and c^Cj is the elec- 
tron density operator. The current is defined as the time 
derivative of the polarization operator given by 



(33) 



In the pure case the current is given by whereas in 

the present case with disorder, the current operator given 
OH 

da 



by ^M- has the following form 



J = ie~ ia / N ^(t + Q2(-l) J )c) +1 c 3 + h.c. (34) 
3 

When the flux a is equal to or n, the system is time- 
reversal symmetric, and there is no persistent current. 
In such cases, we can consider a change of the electric 
polarization by an adiabatic change of parameters Qi 
and Q2. In the linear response theory, the change of the 
electric polarization is given by 



with 



6P = F 1 AQ 2 -F 2 AQ 1 , 



l ^,<*0|J|*m><*m|^|*0> 



Fi = --T( 



<#o|J|# ro )(^ ro |^|^o) 



(35) 



c.c.), (36) 



z.c), (37) 



where L = Na is the circumference of the ring, and a is 
the lattice spacing, and 



OH 
dQi 



E(- 1 ) J 4 c i 



(38) 



dH 1 



dQ 2 



J2(-l) j (e ia/N 4c j+1 + e-^ N c] +1 c09) 



We note that the change of the polarization, r5P, is in gen- 
eral dependent on the path in the parameter (Q) space. 
The |^o) m (|36|) and (|37|) is the many-body ground state, 
and |\£ m ) denote the excited states. Eq and E m are the 
energy for the ground and the excited states, respectively. 
In the pure case when the chemical potential lies in the 
gap, the system is a band insulator. The magnitude of 
the energy gap is given by the magnitude of the CDW 
order parameter (Qi,^)- In the presence of the disor- 
der, the gap closes when the magnitude of the disorder 
potential becomes the order of 0(y/Q\ + Q^)- The sys- 
tem remains insulating because each state is localized, 
i.e., Anderson localization. 

Although we have defined ([35]) , ([37|) and ([36]) only for 
a = and a = 7r, we henceforth extend these formulae to 
general a, which makes the topological properties of the 
charge pumping manifest. We note that except for a = 
and 7r, SP(a) does not mean a change of polarization. 
With this extension to arbitrary a, let us succinctly ab- 
breviate (Qi,Q2joO as Q = (Qi,Q2,Qs) and define the 
gauge potential as 



*(*o|^|*o). 
dQ 



(40) 



This gauge field is so defined that the corresponding field 
strength F = V x A has the components given in ([37]) 
and ([36]) . Furthermore, when the parameters (Qi,Q2) 
are changed along a cycle, the pumped charge can be 
written as 

SP(a) = j dQxa- F(Q) (41) 

where S is a loop on the Q1-Q2 plane. Using the Stoke's 
theorem, ([4T]) becomes 



SP(a) = J d 2 QD(Q) 



(42) 



where D(Q) = d\F\ + $2^2 defined as the distribution 
function of the polarization. This pumped charge is not 
necessarily an integer. 

Nevertheless, in a strongly disordered system it be- 
comes an integer, with an exponentially small correction 
of the order of e~ L ^. To see this, we note that the 
wavefunctions in a strongly localized case are almost in- 
tact with the twisted boundary condition a within an 
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accuracy of e L ^. Hence AP(a) is equal to an average 
AP(a) over a, namely 

AP= [ 27T —AP(a)= [ da-F (43) 

JO 27r JdT 

where dT is the torus surface, and da is the surface el- 
ement. This average AP represents the total flux over 
a torus surface, and it is an integer, which follows from 
topological properties of the gauge field A. As the field 
strength is defined as F = V x A, it satisfies V • F = if 
there is no singularity, i.e. there is no energy crossing. 
Thus when there are no energy degeneracies inside the 
torus V x [0,27r], the averaged pumped charge AP be- 
comes zero due to the Gauss theorem. If there is a energy 
crossing, the gauge field A has a U(l) monopole there, 
and the averaged charge AP becomes a total strength of 
the monopoles inside the torus, again by using the Gauss 
theorem. As the monopole strength is quantized to be 
an integer, AP is always an integer. This discussion for 
([43|) corresponds to that for (2. 23) in Ref. 9, in which 
they claimed that the charge transport is quantized in 
the presence of gap which does not vanish in the ther- 
modynamic limit. On the other hand, the aim of this 
paper targets at another respect: we are interested in 
the strong disorder case, in which the energy gap closes 
in the thermodynamic limit and all wave functions are lo- 
calized characterized by the localization length £ defined 
by £ _1 = J dx\^(x)\ 4z , where ip is the single particle wave 
function. Nevertheless, for finite-sized systems gap al- 
ways exists which is the order of the band width divided 
by the size of the system, which validates the above topo- 
logical argument even in the strongly disordered cases. 
In such disordered systems, the field strength F is highly 
anisotropic around the monopoles, which is closely re- 
lated with a "resonance tunnelling" picture in section 
III.B as we explain in the next paragraph. 

The adiabatic charge transfer is based on the occur- 
rence of the singularity in the parameter space Q, which 
happens when LUMO and HOMO are degenerate. By 
making a close contour enclosing the singularity, the wave 
function shifts as shown in Fig.[lOl Since the disorder po- 
tential is random, we expect that the proliferation of the 
singularities may occur. Furthermore, there could be sin- 
gularity string rather than point in the strong disorder 
limit. The reason is given as the following. In the strong 
disorder limit, all wave functions are localized. There 
could be a situation that the wave function of the LUMO 
and the HOMO are well separated so that the overlap in- 
tegral of two wavefunctions are exponentially small as the 
parameters change. Thus even when the change of pa- 
rameters in numerics or in experiments is very slow, its 
timescale might be still shorter than the inverse of the 
(exponentially small) overlap integral of the LUMO and 
HOMO. Therefore, within this variation of parameters, 
the LUMO and HOMO can be regarded as degenerate. 
Once they are degenerate at a certain point in the pa- 



rameter space, they will keep the near-degeneracy as the 
parameters change until their wave functions have appre- 
ciable overlap and open a gap. As long as they remain 
nearly degenerate as the parameters change, the singu- 
larity string is created in the parameter space. More- 
over, because the wave functions are highly localized on 
the singularity string, they are insensitive to the bound- 
ary condition, in other words, a-independent. Therefore, 
they have the sheet-like structure extending along the 
a-direction in the parameter space. Even though it is 
sheet-like, we call it a string referring to its projection on 
the (Qi,Q2) plane. From this observation, we conclude 
that the singularity in the parameter space projected on 
Qi and Q2 subspace has two kinds, points and strings, 
in the strong disorder limit. The former contribute to 
the charge transfer, while the latter do not. The charge 
transfer for the former is explained within the resonance 
tunnelling Section . 

Numerical Result 

In this section, we solve numerically the eigenvalue 
problem for slowly varying parameters (Qi, Q2) and cal- 
culate the contour integral in flS]) . We introduce the 
small segment AQ for the contour integral, and demand 
that Fi,2<2 are continuous. Therefore, the size |AQ| de- 
termines the energy-scale and time-scale of our simula- 
tion. In actual calculation, we take |AQ| = 10 -6 t n .n. 
as the finest one. Then the discussion below applies 
when we observe the system within the time-scale T ~ 
h/\AQ\ ~ 10 6 ft/t n . n .. Also we average over a to see the 
quantization of the charge transfer in (|43|) . 

We solve (|3T|) by a direct numerical calculation. Our 
goal is to compute the distribution function from the 
(|42|) . which is a-dependent. We will also demonstrate 
the quantization of the charge transfer by averaging over 
a by (1431) . 

We consider a ring with N = 50 and take Vi/t n , n , to 
distributes uniformly in [—7.0, 7.0]. To obtain the distri- 
bution function, we divide the (Qi,Q2) plane into small 
pieces of grids and calculate (|42|) for each grid. Let us 
remind the readers that the distribution function for the 
pure case is the S— function at the origin[6]. In that case, 
we do not need to integrate over a to obtain the quan- 
tization of the charge transfer, because the singularity 
lives in the (Qi, Q2, k) space, and the integration over k 
is already included implicitly in (|42|) . 

Our result can be summarized by the cartoon in 
Fig. [JJ] First, we found that AP = 1.0 if we take the 
loop A which is the outer dash one (|Q| ~ 10) and av- 
erage over a. It corresponds to the weak disorder case 
where the strength of the disorder is much smaller than 
|Q|, so the potential strength is not large enough to 
close the gap. According to Niu and Thouless[9], the 
polarization is quantized, and we reproduce their result 
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here. Secondly, the singularity strings appear in the re- 
gion \Q\ < 4.5. In this region, the disorder potential is 
strong enough to close the gap. As we pointed out in the 
section , there are two kinds of the singularities, points 
and strings. In the "stringful" region, it is rather dif- 
ficult to locate the singularity point. It is because the 
operation of the adiabatic process changes when making 
the contour across a string. In Fig. [12j we illustrate two 
different adiabatic processes regarding making contours 
enclosing the singularity points and across the singularity 
string. Fig. fT2l a) shows the usual definition of the adia- 
batic process. In the presence of the gap, the adiabatic 
process has to follow the lower energy state, which is in- 
dicated by the square dots in Fig. H^f a). Fig. [T2lf b) shows 
the adiabatic process of the contour across the string. In 
this case, because the wave functions of the HOMO and 
the LUMO do not overlap, the adiabatic process needs 
to keep ramping up to the higher energy state at the 
first degenerate point and ramp down to the lower en- 
ergy state at a second degenerate point, which completes 
a close loop. Note that the energy of the HOMO is larger 
than that of the LUMO after the first degenerate point 
and becomes smaller again after the second degenerate 
point. In other words, the system stays in the same state, 
and thus the charge does not transfer at all after a close 
contour. Therefore, even if we make a close loop which 
encloses a whole singularity string, the charge does not 
transfer. It should be noted that even though the sys- 
tem is at a higher energy state, it takes ~ e L ^ to relax, 
where L is the size of the system, and £ is the typical 
localization length which is around 3 sites in our case. 

In Fig. [T2lf c). we demonstrate a real process across a 
singularity string indicated by the loop C in Fig.[TTl Usu- 
ally, we refer singularity string to the degenerate mani- 
fold between the HOMO and LUMO. In a real process, we 
may have to consider the degeneracy between the HOMO 
and the 2 nd LUMO. For example, after the degenerate 
point between LUMO and HOMO, the descending order 
in energy is 2 nd LUMO, HOMO, and the LUMO. There is 
the possibility, for example the boundary (3 in Fig. [T2l c). 
that the energy of 2 nd LUMO and that of HOMO can 
be the same along the contour before meeting the next 
degenerate point between HOMO and LUMO. If it hap- 
pens, the system still stays in the original HOMO state 
by ramping further up to the higher energy state. After 
that, the descending order in energy becomes HOMO, 2 nd 
LUMO, and LUMO. Fig. fT2Tc) shows the complete pro- 
cess and the system has to go back to its original state af- 
ter completing a close loop. As a result, HOMO, LUMO, 
and other states can be viewed as the Riemann sheets. 
Singularity strings are like the branch cuts which con- 
nect different Riemann sheets. Contours on the (Qi, Q2) 
plane are very similar to the one in the complex plane. 

As the string just corresponds to the degenerate states 
which does not penetrate through the whole sample, 
the only topological object which contributes to the po- 



larization is the isolated singularity point We found 
one at (Q 1 ,Q 2 ,at) = (-0.3885,-4.80415,0). If we cal- 
culate ([43]) around that point, for example loop B in 
Fig- HH we obtain AP = —1.0. We identify the sin- 
gularity point which carries AP = —1.0 with the an- 
timonopole in contract to the monopole which carries 
AP = 1.0. Considering together with the total flux 
calculated from the loop A in Fig. [TT] to be 1.0, there 
must be at least two monopoles and one antimonopole 
in our system, which demonstrates the proliferation of 
the monopole- antimonopole pairs in the strong disorder 
limit. Unfortunately, we are not able to locate monopoles 
exactly but obtain rough position, since they merge in 
the string- ful soup. The integral in (jHj) converges very 
slowly as the counters are close to the singularity string 
(or points). The integral step is equivalent to the rate 
of the adiabatic process, which has to be less than the 
typical energy spacing 2(£ + v* nax )/A/\ When approaching 
to the singularity, the energy gap becomes very small, so 
the process rate has to be smaller to have convergent re- 
sult. Therefore, it is very difficult to locate the monopole 
position if they are in the string-ful soup. 

a 2 




A 



FIG. 11: A cartoon picture of the distribution of singulari- 
ties at a = 0. The singularity has two kinds: strings and 
points. Calculating (|41[) from the outer big dash loop gives 
the polarization +1.0 and -1.0 from the loop B. If we stay 
far away from the origin, one would think that there is only 
one monopole with AP = 1.0 at the center like the pure case. 
However, when the disorder strength is strong enough, there 
are some subtle structure inside the big monopole. There are 
the monopole- antimonopole pairs and the singularity strings. 

The appearance of the antimonopoles in the strong dis- 
order limit is very appealing. It contributes to the po- 
larization in the opposite way to that a monopole does. 
The polarization comes from the charge transportation 
through the whole sample from one end to the other. 
Therefore, (|42|) for the antimonopole must be highly a- 
dependent shown in the Fig. fl3lfa). The peak value of 
AP(a) happens at a = 0. After integrating over a by 
(|43|) . we obtained AP = —1.0. This result indicates that 
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a = 0. Fig. [T3lf e) shows the energy gap along the valley. 
The slope of the energy gap function along the valley is 
roughly 10 times smaller than other directions. In other 
words, the antimonopole looks like a football more than 
a round ball. This anisotropy can be understood in terms 
of the resonance tunnelling as we see in the next section. 
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FIG. 12: The adiabatic process in the gap less case is different 
from the gap case, (a) shows the process with a gap. (b) 
shows the process without a gap. (c) a real process from loop 
C in the Fig. [II] 



the antimonopole may be located at a = 0. In Fig fl3lf b). 
we found that the energy gap closes at a = supporting 
this observation. In Fig. fl3lf c) and fl3lf d). we show the 
energy gap at a = along Qi and Q2 direction, respec- 
tively. In fact, the structure of the antimonopole is highly 
anisotropic. The energy gap has a valley-like structure at 
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FIG. 13: Some properties of antimonopole. (a) the a- 
dependence of (|42p around the antimonopole. (b), (c), and (d) 
the energy gap functions of the antimonopole along the a, Qi, 
and Q2 direction, (d) the gap function along the valley. The 
antimonopole is not spherically symmetric. The gap function 
at a = has a deep valley because the slopes along the val- 
ley is 10 times smaller than other directions, (f) The wave 
functions of the HOMO and LUMO at the antimonopole. In 
Figs, (b)-(e) the axes are in the unit of t n .n.- 

Another quantity to measure the extendedness of a 
state is the Thouless number defined by 



M T = 



E(a = tt) - E(a = 0) 
SE 



(44) 



where E(a) is the energy of the state and SE is the typ- 
ical energy spacing. We plot the Thouless number of the 
HOMO in Fig. [TH We found a ridge distribution which 
coincides with the valley in the energy gap function. The 
antimonopole locates at the ridge but not necessarily the 
highest one. The Thouless number at the ridge is one 
order of magnitude larger than the other points in the 
vicinity in the parameter space. It suggests that the 
high boundary condition sensitivity must dominate the 



14 



physics at the antimonopole. It leads to the resonance 
tunnelling in the Anderson insulator. [33] 




FIG. 14: The Thouless number in the vicinity of the anti- 
monopole. The Thouless number has a ridge structure which 
exactly coincides with the valley in the gap function. 



Resonance Tunnelling 

We now compare the numerical results with the picture 
of resonance tunnelling discussed in Sec. . In the reso- 
nance tunnelling scenario, we can deal with the periodic 
system by identifying xs with x\. It leads to k$ = k±. In 
such a periodic system, as is apart from an open system, 
Ep is no longer a tunable parameter, but corresponds to 
discrete energy levels of the system. Thus, S^, c^, and 
fa are smooth functions of Qi, Q2 and an energy E. Be- 
cause the periodic system is threaded by a flux of a<I>o, 
we require that 

--(£)• <«» 

It yields det(0 — e l0i ) = 0, i.e. Re#n = cos a. Because 
On is a function of Qi, Q2 and E, this defines a (discrete) 
energy level of the system as a function of Qi and Q2. 

We study the monopoles in the framework of resonance 
tunnelling. As the monopoles appear at band-crossing 
points, we have to look for degenerate solutions of ([45]) . 
Namely, we look for cases with two (degenerate) states 
for fixed Qi, Q2 and E. This happens if and only if the 
matrix is equal to e za . This condition reduces to 

Si = S 2 , (46) 

uj^2 k(x)dx = (2n + 1)tt, (47) 

uo' = 2 / k(x)dx + 2 k(x)dx = (2n + 1)tt, 



(48) 

a = -{n + n 1 + 1)tt, (49) 

where n and n' are integers. Because Si, S2, uo and 
uo' are functions of Qi, Q2 and E, the three conditions 
(fl6]h (jUj) and ([15]) determine a set of isolated points in 
the Q1-Q2-E space. Combining with ([49]) , we get iso- 
lated points in the Q^-Q^-a space, corresponding to the 
monopoles/antimonopoles [38j |. Except for the vicinity 
of the monopoles/antimonopoles, the field F is almost 
a-independent, because the system is almost intact with 
the change of a. Thus the distribution of the field F 
is as shown in Fig. [T5l a). In contrast with a small- 
sized model £ ~ L in ref. 38, in our disordered model 
(£ <C L) the "near-field" region where the flux density is 
a-dependent is limited only in the very close vicinity of 
the monopoles/antimonopoles. 

Equations (jUj) and ([35]) are the Bohr quantization con- 
ditions, meaning that the there are two localized states 
around X2 and £3 (= x±) : having the same energy. In gen- 
eral, when there are two localized states with the same 
energy, displaced each other by a high potential peak, 
there is a small tunnelling matrix element between them. 
It gives a small energy splitting between bonding and 
anti-bonding states. As a result, the degeneracy is lifted. 
Nevertheless, when the equality ([46]) holds, it guarantees 
that the hybridization between them cancels exactly, and 
the degeneracy is not lifted even after one takes the tun- 
nelling into account. Thus the two localized states are 
in resonance with each other. Equation ([49]) means that 
such exact degeneracy occurs only when a is equal to 7r 
or 0, confirming our numerical result. 

If we change Qi and Q2 to make Si 7^ S2 with the 
conditions (|47j) and (|48|) preserved, there occurs a small 
splitting (~ 0(e~ 2Si )) to the otherwise degenerate states, 
due to a small unbalance between Si and S2 . These con- 
ditions (|I7j) and (|38|) specify one direction in the Q1-Q2 
plane from the monopole. Along this direction the gap 
remains nonzero but very small, thereby the field F be- 
comes large as schematically shown in Fig. [15] (b). In 
our numerical results, this direction corresponds to the 
valley direction in the energy gap, and the ridge in the 
Thouless number in Fig. [E] Thus the anisotropy of the 
energy gap in the Q1-Q2 plane is of the order e 2Si ~ e L ^. 
Because in our numerical calculation L = 50, £ ~ 7, the 
anisotropy of the monopole is e 2Si ~ e L ^ ~ e 7 ~ 10 3 , 
which should be compared with the numerical value ~ 10 
in the previous subsection. Their difference can be at- 
tributed to finiteness of the mesh size in the numerical 
calculation; the mesh may not be fine enough to repro- 
duce the anisotropy ~ 10 3 . We also note that the strings 
observed in the numerical calculation can be regarded as 
the curve of uo = (2n + l)7r, uj' = (2n' + l)7r, with Si and 
S2 are unrestricted. 
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o=(2n+1)n 
co =(2n' +1}jc 



FIG. 15: Distribution of the vector field F around the 
monopole. (a) the field F is mostly uniform in the a-direction 
because ^ < L, except for the vicinity of the monopole. (b) 
the field is anisotropic in the Q1-Q2 plane. 



Discussion 

In Fig.([T3f). we show the wave functions of the HOMO 
and LUMO, which supports the resonance tunnelling 
scheme. At the resonance, the energy of HOMO and 
LUMO are degenerate. If HOMO sits at x u LUMO can 
be viewed as the resonance state at the potential valley 
#2- Therefore, by means of LUMO, HOMO can tunnel 
through the whole sample. We must emphasize that this 
mechanism is topological as that in the pure case. 

The adiabatic rate depends on the energy gap. At 
the resonance, the wave functions of the HOMO and 
the LUMO do not overlap. By tuning the parameters, 
the overlapping integral between these two states are not 
zero, and thus the gap opens up. Again, the rates that 
the gap opens with the parameters are not necessarily the 
same. Therefore, the monopole (antimonopole) struc- 
ture is not necessarily isotropic. If fact, it can be highly 
anisotropic as shown in our calculation. Then, the adi- 
abatic rate is restricted by the minimal gap along the 
contour. It can be roughly estimated as e~ L ^ R , where 
£r is the localization length at the resonance. 

The current results for the periodic/twisted boundary 
condition is consistent with those for the open boundary 
condition. The topological structure in the open bound- 
ary condition, namely the vortices, can be regarded as a 
projection of the monopoles to a 2-dimensional plane. 
The anisotropic flux distribution and the anisotropic 
phase gradient are two facets of the resonance tunnelling. 
While the results for the open boundary condition is 
more applicable to the real materials and design, the non- 
trivial generation of the monopole-anti-monopole shown 
in the periodic/twisted boundary condition case illus- 
trates the deep and whole new topological structure in 



the strong disordered systems. 



CONCLUSION 

In this paper, we have studied the charge pumping and 
dielectric response in disordered insulator for both the 
open boundary condition and periodic/twisted boundary 
condition. As for the open system coupled to the leads, 
we have found the quantum-mechanically enhanced di- 
electric response in nanoscopic/mesoscopic disordered in- 
sulators, which give a guide for fast and low-dissipation 
ferroelectric thin films of the FeRAM. The phase of the 
reflection coefficient r is a key parameter for the pump- 
ing, and has the rich structure in the parameter plane. 
Topological nature of the insulator dictates the phase 
winding of r, around the vortex where r = and \t\ = 1. 
In a pure insulator, this corresponds to the gapless point. 
In the disordered case, it corresponds to the resonance 
tunnelling through the sample, whose position varies 
with the chemical potential. 

With the periodic closed system, we can consider 
the three dimensional parameter space Q = (Qi,Q2,<^) 
where Qi is the site energy alternation, Q2 is the bond 
dimerization, and a is the phase twist for the boundary 
condition. In this three dimensional space, one can dis- 
cuss the monopoles and the their associated gauge field 
distribution, which are again interpreted in terms of the 
resonant tunnelling. 

The existence of the resonant tunnelling in the param- 
eter space is guaranteed by the topological properties of 
the charge pumping, and only two parameters such as Qi 
and Q2 need to be tuned to realize it. These parameters 
can be controlled experimentally by the external electric 
field and pressure. Therefore the present theory offers a 
way to design the enhanced dielectric response in realistic 
systems by controlling the disorder. 
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